Influence of mechanical properties on mixing and segregation of granular material 
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Processes of mixing and segregation in a packed bed of granular material stirred by a period- 
ically moving rectangular bar are simulated by the discrete element method (DEM). Influence of 
mechanical properties of the particle material (dynamic friction coefficient, density, elastic modulus) 
on the segregation intensity is researched. Dynamic friction coefficient was found to have the most 
noticeable influence to the segregation intensity. 



I. INTRODUCTION 

Granular materials of various forms are ubiquitous in 
nature and technology. They exhibit a rich variety of 
phenomena and, depending on circumstances, can have 
properties similar to either solids, liquids or gases [ij, or 
behave in a completely different manner. Size segregation 
is a characteristic feature of the moving granular media 
and it is important, as well as mixing, in industrial 
handling of granulated materials. A number of physical 
mechanisms were suggested to explain the segregation 
process Q. Existing theories of continuous mechanics 
or statistical physics have only limited applications for 
description of granular media, and a unified theory en- 
coinpassing all the granular phenomena is still missing 
d, Q . An ability to describe and predict the conduct of 
granulated materials is important for technological ap- 
plications in industry where bulk materials are routinely 
handled. On the other hand, granular media is an in- 
teresting object for theoretical studies because of its in- 
tricate phenomenology despite ostensibly simple purely 
mechanical nature. Numerical simulations can provide a 
useful insight into the origins of various aspects of this 
behaviour. 

Research data about the influence of the mechanical 
properties of granular matter to the segregation process 
is still scarce. Besides, simplified models are commonly 
used: binary mixtures of particles are researched, binary 
collisions between the particles are assumed. Experi- 
ments involving the particles with different mechanical 
properties are difficult to implement, but numerical sim- 
ulations with different parameter sets can be performed 
easily, even though they are lengthy and require much 
computing power. 

The objective of the present paper is to investigate the 
mixing process inside granular material and the depen- 
dency of mixing on the properties of the particle material. 
Most research performed in segregation of granular ma- 
terials dealt with the situation when the researched ma- 
terial is put to motion by a moving bottom wall Q . We 



modelled a different setup where the material is stirred 
by a rectangular bar moving forward and backward in 
the horizontal direction and buried inside the material, 
which is a model of a simplified case of real equipment, 
such as forward or backward moving grates widely used 
in industry. Mixing of granular material is sustained by 
periodic motion of a rectangular bar inside a rectangular 
container. We observe the segregation and stratification 
of particles by their sizes and intensity of mixing. 



II. SIMULATION MODEL 

A number of models have been applied for describ- 
ing granular media with a varying degree of success . 
The discrete element method (DEM) proved to be the 
most accurate, besides, it is simple to implement. In 
this method, the motions and collisions of each separate 
particle are tracked using the equations of classical me- 
chanics. The equations of motion are integrated with a 
constant step size considerably shorter than the duration 
of the collision. Particle motion obeys the usual Newton 
equations 
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where is the i-th particle position, is its velocity, 
Fi = ^ Fij + TTiig is the total force acting upon the i- 

th particle, F^j is the force acting upon the i-th particle 
arising from its collision with the j-th particle, mi is the 
particle mass, g is the free-fall acceleration vector, is 
the angular velocity, = ddj x F^ is the torque acting 
upon the particle, d^j is the collision point vector defined 
below, li is the moment of inertia which is scalar for 
spherical particles, N is the total number of particles. 
The orientations of the particles are not updated because 
this has no sense for spherical shape. 

The scheme of collision between two particles is shown 
in Fig. [TJ Particle deformations during collisions are ap- 
proximated by a partial overlap of their geometric shapes, 
with the overlap depth hij . The position of colliding par- 
ticles i and j in space and their motion is defined by the 
radii- vectors Xi, Xj of their centres Oi, Oj, linear veloc- 
ities Vi, Vj and angular velocities w^, wj. The relative 



position of the colliding particles is 



The 



point of contact Cij is located in the middle of the over- 
lap area and its radius- vector is y.cij ■ The unit vector riy 
is normal to the contact surface and directed toward the 
particle i, n^j — — n^i. The vectors ddj and dcji specify 
the position of the contact point dj with respect to the 
centres of the particles i and j: ddj = — x^. Relative 
velocity of the particles at the contact point is defined as 
= - Vcji, = -Vji, where Vc^j = Vi-t-w; x d^j, 
^cji — Vj + Wj X d^ji are the velocities of the parti- 
cles i and j at the contact point. Normal and tangen- 
tial components of the relative velocity are expressed as 
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respectively. 



The contact forces act at the contact point. The collid- 
ing particles can slip with respect to each other in the 
tangential direction by the distance St^ij = |/ vt,ij(t) dt| 
which is referred to as the integrated slip. As the parti- 
cles slip, the value of 5t^ij increases until the tangential 
force exceeds the force caused by static friction, and af- 
ter that, the tangential force component is defined by the 
dynamic friction. 

There is still no definite agreement what force model 
best describes the real interaction of the colliding par- 
ticles @, 0]- Two most popular models of the normal 
forces arising between the colliding particles i and j are 
the Hook's law and the Hertz's law. These laws are ex- 
pressed respectively as 
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It can be expected that the effect of a particular force 
model will be "averaged out" over the vast number of 
collisions during the process and the differences between 
these models will not have a significant effect on the 
statistical parameters of the granular matter. In order 
to verify this, we performed the simulations using both 
models and compared the results that will be presented 
below. 

The tangential component of the interaction force is 
expressed as 




FIG. 1: Geometrical scheme of the elastic interaction 
between the colliding particles 



the shear modulus. The simulations described below used 
the time-driven DEM implementation . The equa- 
tions of motion ((!]) were integrated using the G^'^-order 
Gear predictor-corrector scheme jH. In order to reduce 
the amount of resources required for collision detection, 
the simulation area was partitioned into rectangular cells, 
with the length of the cell edge equal to the largest par- 
ticle diameter. During each time iteration, initially all 
the particles were assigned to the corresponding cell ac- 
cording to their current coordinates (Verlet lists). It is 
evident that only the particles belonging to the same cell 
or the nearest neghbouring cells can collide; thereby the 
number of pairwise collision tests can be reduced from 
0{N^) to almost 0{N). 
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E is the elastic modulus of the 
r, ■ To 



particle material, v is the Poisson's ratio. 



iriij — — are the normalized radius and normal- 
mi rrij 

ized mass, hij is the overlap depth of the particle shapes 
during collisions jB], ny and Uj are the normal and tan- 
gential vectors at the contact point, 7n and "ft are the 
normal and tangential (shear) dissipation coefficients, fi 
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is the dynamic friction coefficient, fct = — ■ , G is 
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III. CHARACTERIZATION OF MIXING AND 
SEGREGATION 

The approach to characterization of mixing of particles 
in granular media stems from the theory of turbulence in 
fluid dynamics fTo'|. It is assumed that the velocity of 
each particle can be split into its mean value and 
the fluctuating part v^: 



v:(t)=v,-V,(i). 
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The mean velocity is defined by averaging the velocities 
of particles surrounding the given i-th particle within a 
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certain preset volume: 
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Here Rv is the radius of the selected volume surrounding 
the z-th particle over which the averaging is carried out, 
Xi is the position vector of the i-th particle, and a logical 
function S is introduced for brevity: 



J (condition) 



1, if condition is satisfied 

0, otherwise. ^ ' 



The fluctuating velocity shows how quickly the i-th par- 
ticle changes its position with respect to the surrounding 
particles. By averaging v'^^ and v'^ over the simula- 
tion volume, the time-dependent mixing parameters are 
introduced: 
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where the coordinate axes are denoted as x = xi, y = X2, 
z = X3. Averaging over the simulation time period [ti, ^2] 
yields the mixing parameters during this simulation pe- 
riod: 
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The segregation is estimated quantitatively as the 
change of vertical positions of the particles during the 
simulation period depending on particle sizes. Particles 
whose radii are within the upper 20% of the total inter- 
val of the particle radii are considered "large" , and those 
with the radii within the lower 20% of the total interval 
are considered "small", i.e., the interval of radii of small 
particles is [0.005 m, 0.007 m], and that of the large par- 
ticles [0.013 m, 0.015 m]. The average vertical position 
of particles whose radii are in the interval ri < < r2 is 
defined as 
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The changes of the average vertical position (y^) of 
"small" and "large" particles are compared to the change 
of the average vertical position of all the particles (includ- 
ing large and small ones). However, as the stirring bar 
moves forward and backwards, the height of the volume 
occupied by the particles ( "packed bed" ) changes. There- 
fore, it is reasonable to introduce the average height of 
this volume HpB{t) and to normalize the vertical posi- 
tions of the particles to this height, in order to discrimi- 
nate the vertical motion of individual particles from the 
motion of the entire packed bed as a whole: 



HpB{t) ' 
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IV. SIMULATION RESULTS AND DISCUSSION 

The simulated system consists of a rectangular box 
whose size is 0.2 x 0.4 x 0.2 m^ containing N — 1600 
spherical particles of various sizes with random distribu- 
tion of particle depths for particles of different sizes, with 
radii ranging from 0.005 m to 0.015 m. At the lower left 
corner of the box, a rectangular moving bar is located 
whose width (in the direction of the z axis) is equal to 
the width of the box (Fig. [2]). The histogram of particle 
size distributions is shown in Fig. O During the process, 
the bar moves from its leftmost position to the right along 
the X axis over the distance of 0.1 m at a constant veloc- 
ity for the time period T, then retracts back to its initial 
position during the same time period T, and then the 
process repeats periodically during the entire simulation 
time. The time of motion of the bar from the leftmost 
to the rightmost position was usually T = 10 sec, if not 
specified otherwise; the time of backward motion is equal 
to the time of forward motion. The bar motion results 
in stirring of the particles and their redistribution in the 
volume of the box. 

The simulations were performed for particles that had 
initially the same initial positions, linear and angular ve- 
locities for each run. The "standard" mechanical param- 
eters of the particle material are listed in Table [H only 
the value of the parameter whose influence was to be de- 
termined was set different for the corresponding run. The 
set of equations H]) was then solved for a certain period 
of the "simulated" time, yielding the positions and veloc- 
ities of each particle at the time moments separated by a 
certain time step. The duration of the simulation period 
usually varied from 100 s to 400 s long. From the result- 
ing data, the mixing and segregation parameters ([5]) - 
(fTTj) were calculated. 

Segregation of particles can be seen by comparing par- 
ticle positions at different time moments in Fig [51 Ev- 
idently, smaller particles tend to sink downwards and 
accumulate near the bottom as the process progresses. 
Quantitatively the mixing can be estimated from the evo- 
lution of the mixing parameters n(t) (Eq. as shown 
in Fig. m 
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FIG. 3: Distribution of particle radii in the simulated 
system. 



TABLE I: Standard mechanical properties of the par- 
ticle material 



Parameter 


Notation 


Value 


Minimal particle radius 




0.005 m 


Maximum particle radius 


^max 


0.0f5 m 


Density 


Q 


700 kg/m^ 


Elastic modulus 


E 


fO^ Pa 


Poisson modulus 


a 


0.2 


Normal dissipation coefficient 


7n 


too s^^ 


Shear dissipation coefficient 


7t 


too s^^ 


Dynamic friction coefficient 


M 


0.5 


Shear modulus 


G 


3 X fO*' Pa 


Period of the stirring bar motion 


T + T 


fO+fO s 


Free fall acceleration module 


|g| 


fO m/s^ 



FIG. 2: Evolution of particle positions during the mix- 
ing process (top to bottom): initial positions ait — sec, 
positions at f = f sec, and at t = 400 sec. Particles of 
different diameters are shown in different colors. 



Segregation is demonstrated by the change of the nor- 
malized average vertical positions {y{t)) (Eq. fTTj) of small 
and large particles (Fig. [5]). It shows evidently how the 
small and large particles separate vertically during the 
mixing process. The jagged shape of the curve corre- 
sponds to upward and downward motion of the particles 
following the motion of the bar. As the process pro- 
gresses, the particles tend to approach a certain nearly 
stationary state where the fluctuations of their positions 
are influenced only by the motion of the stirring bar. 
These fluctuations are more pronounced for small parti- 
cles accumulating near the bottom, because the available 
volume changes as the bar moves. The particles accumu- 
lated at the top of the packed bed are located mainly 
above the bar and are therefore less influenced by its 
motion. 

The time-averaged mixing coefficients Pb, PB,x,y,z 
(Eq. [S]) are calculated from the time-dependent mixing 
coefficients Tl{t), Tlx.y,z{t) by averaging over a certain 
selected time period [ti,t2]. As seen from the tempo- 
ral evolution of the vertical particle positions, in case 
of the random initial distribution of particle positions 
by their size, mixing is most intense in the beginning 
of the process, while particles are not yet segregated by 
their sizes. Therefore, the values of the parameters Pb, 
PB,x,y,z depend on the selection of the averaging period 
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FIG. 4: Evolution of the mixing parameters U for the piG. 7: Evolution of vertical positions of smaU (curves 

values of the dynamic friction coefficient = 0.2 and 4)^ j^rge (curves 2, 5) and aU (curves 3, 6) particles 

^^ ~ ^-^^ with the densities of q = 500 kg/m'^ (curves 1, 2, 3) and 

Q = 900 kg/m^ (curves 4, 5, 6). 

























fi = 0.2 - small 

H = 0.2 - large 












fi = 1 - small 

M- = 1 - large 























50 100 150 200 250 300 350 400 



FIG. 5: Evolution of the average vertical positions of 
small and large particles for the values of the dynamic 
friction coefficient fi = 0.2 and fi = 1.0. 
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FIG. 6: Dependence of the averaged mixing parameters 
Pb on the dynamic friction coefficient p. 



[ii,t2]- It is possible that some parameters of the mate- 
rial or the mixing procedure have more influence during 
the initial stages of the process, while other parameters 
are more significant as the material is approaching the 
steady state. 

As seen from Figs. |3]-[5l higher values of the dynamic 
friction coefficient causes more intense mixing. This re- 
sult is similar to the case of hydrodynamics where mixing 
is more intense in more viscous fluids. The dependence 
of the overal mixing intensity on the dynamic friction 
coefficient is nearly linear in the researched range of its 
values (Fig.[n]). 



In order to identify the influence of the material den- 
sity of the particle material, simulations were performed 
for particles with the densities g = 500 kg/m'^ and 
g — 900 kg/m"^. The changes of the vertical particle po- 
sitions for these different values of the density is shown 
in Fig. [71 Evidently, the overal effect of the material den- 
sity on the segregation process is insubstantial. However, 
it can be noticed that at the later stages of the process 
the larger particles with lower density accumulate slightly 
higher in the packed bed than in the case of denser par- 
ticles. This results in looser packing of the particles in 
the top layers of the bed. 

Fig. [5] demonstrates the evolution of the vertical posi- 
tions of particles interacting according to the Hook's law 
(Eq. [2]) compared to the case of interaction according to 
the Hertz law (Eq. As noted above, the difference 
between the two force models does not manifest itself 
signiflcantly in the mixing process. The simulations us- 
ing the Hertz normal force model were performed with 
the values of the dynamic friction coefficient fi ~ 0.5 and 
/i = 0.8. Table [IT] lists the values of the mixing coefficient 
Pb obtained for these different values of the dynamic 
friction coefficient for both Hook and Hertz models. The 
values of Pb were calculated by averaging U{t) over the 
time interval [1,100] s, when the mixing process is still 
more intense than in the later stages. The difference be- 
tween the Hook and Hertz models yields the difference in 
the mixing coefficient of approximately 10%. It can be 
noticed that the difference of this coefficient for different 
values of the dynamic friction coefficient is larger in case 
of the Hertz model; however, more simulations with dif- 
ferent parameter sets are required in order to determine 
the sensitivity of the Hertz model to the mechanical pa- 
rameters. 

Dependence of the mixing process on the elastic mod- 
ulus was also investigated in the range E — (0.6 -f- 1.4) x 
10"^ Pa. The values of the mixing parameters Pb, PB,x,y,z 
decrease slightly as the values of E increase, however, this 
change was insignificant, within 2%. On the other hand, 
the range of the values of E was quite small and the ef- 
fects of this parameter were probably masked by other 
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FIG. 8: Evolution of normalized average vertical posi- 
tions {Y{t)) of small and large particles interacting ac- 
cording to the Hook's law and the Hertz's law. 



TABLE II: Comparison of the values of the mixing 
coefficients Pb (in cm/s) averaged over the period from 
ti = s to ^2 = 100 s for the cases of interaction accord- 
ing to the Hook's law and the Hertz's law for different 
values of the dynamic friction coefficient: fi = 0.5 and 
fi = 0.8. 



Normal force model 


fi = 0.5 


fi = 0.8 


Hook 


1.3569 


1.6780 


Hertz 


1.2947 


1.7335 



V. CONCLUSIONS 

The mixing parameters introduced by the hydrody- 
namic analogy were applied to characterize the mixing 
and segregation of granular matter. Segregation of par- 
ticles of different sizes was estimated by the change of 
their normalized depths. 

The mixing intensity weakly depends on the elastic co- 
efficient of the particle material in the researched range of 
the values. In order to better understand the influence of 
this parameter, it is necessary to research a wider range 
of its values. On the other hand, the mixing intensity 
depends noticeably on the value of the dynamic friction 
coefficient. In turn, more intense mixing results in faster 
segregation. This confirms the validity of the mixing pa- 
rameters for estimating the segregation. The dependence 
of the mixing intensity on the dynamic friction coefficient 
is linear in the researched range of the values. 

There is a certain difference of the obtained results de- 
pending on whether the Hook's or Hertz's model of the 
contact force between the colliding particles is used; how- 
ever, this difference is not substantial. Further research is 
needed to identify possible influence of some parameters, 
by extending the range of the researched values. 
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